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ABSTRACT 


A predictor correction technique to improve the performance 
of Weapon Control Systems has been developed and found to work 
satisfactorily for various . typically chosen cases. The logical 
organisation of a special purpose Digital Computer to be used vrith 
the Weapon Control System is described. 

This is organised in two parts. In the first part the 
predic'tor correction technique is developed and in the second part 
the logical organisation of a special purpose computer is described. 

Chapter 1 gives the general picture of the fire control 
problem that is solved by the predictor in the Weapon Control System. ; 
In Chapter 2, the possible sources of error in prediction are discussed.' 
In Chapter the error in prediction is analysed and an error correctioi 
technique is deviced. In Chapter 4> the logical organisation of a 
special purpose computer is described. A gate level design of the 
computer is given and its performance compared with an existing 
electro-mechanical computer. 



PRSDIGTCE CCItSaCTICH TSCMI-^US 



CHfi-PEER 1 


IHTRODUCTIQtT 


1.1 Greaeral 

hi weapon control systems dealing with mobile targets, a 
predictor is used to compute the angle at which a projectile is to be 
fired so as to destroy Hie moving target. This predictor is being 
provided with information regarding the present position and rate of 
change of position of the target. The predictor assumes that the 
target, will maintain its present course, that is, it will continue to 
move in the same direction at the same speed as now, and predicts its 
future position and the angle at which a projectile is to be fired so that 
the target and the projectile being fired at this mcsn^nt, will collide 
after certain time Interval. Bais time interval, is called the "time 
of flight of the projectile". 

The future position of the target and the time of flight of the 
shell are interdependant. Hence the predictor in effect solves a set of 
simultaneous equations. Row a simplified weapon control system using a 
predictor will be considered. 

1 .2 General Representation of the System 

A simplified weapon ccntrol system used to control an anti- 
aircraft gun is shown in Fig. 1 . 1 . To simplify the system further, the 
following assumptions are made in the beginnings 

1) 353.6 gtn and Radar interna are mormted on a stabilised 


platfoimi 




B1 ook Dirgram of a Simplified Weapon Control System 
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2) The Gm and the J?adar Antenna are theoretically 
at the same point and plane. 

These assumptions will be taken care off in the following 
manner. The first assumption is taken -care off by applying the correcting 
signals from a Gyro Stabiliser and the second assmption by applying 
lateral and vertical convergence corrections to the angle at T.diich the 
projectile is to be’ fired. It is also assumed that the aircraft will 
continue to fly in the same course. This will be justified in a later 
section. 

'When, an aircraft is detected, the Radar will automatically track 
the aircraft. Kecessary signals for this automatic tracking are generated 
by -Hie Tracker. Ihile tracking' the aircraft automatically, the following 
informations are obtained* 

1 ) Itesent range of the aircraft, R.^. 

2) Eresent bearing of the aircraft, T . 

5) Ecesent elevation of the aircraft, S^, 

4) Rate of change of range of the aircraft, 

5) Hate of change of bearing of the aircraft, 

6) Rate of change of elevation of the aircraft, 

These signals are fed to the predictor. The wind velocity and 
direction of wind are also fed to the predictor to apply necessary wind 
corrections. The muzzle velocity, that is, the velocity of the projectile 
at the time of leaving the Gun, and the ballistics of the S fee 11 which form 
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another set of inputs to the predictor are manually set. The muzzle 
m velocity and ballistics vary for different types of Guns and Shells. 

Erom these inputs the predictor is req^uired to compute the followings 

1 ) -Angle in the horizontal plane to which the Gun 

should be trained before firing a projectile, T 

g 

2) Angle of elevation of the Gun, E 

S 

5) Time of flight of the Shell, t^ to set the 
fuze for the Shell. 

The ftmction of the predictor is to compute E^, T^ and t^. from 
the available inputs. This problem is called as the fire control problem. 

This is represented by a set of simultaneous equations. These equations 
are not given here as they resemble an existing iwapon control system* 

1 .5 -i^ossible Sources of Error 

In this weapon control system, to simplify the function of the 
predictor, many approximations are made. In the set of simultaneous | 

equations used to solve the fire control problem these are incorporated. 

These equations involve certain constants whose values depend upon the ; 

bearing, altitude , speed and range of the aircraft and the wind conditions 
at that tine. 

Though in a real situation they are not constants but for simplifying 
the predictor, they are assumed to be constants. They are set at their 
optimum values which are determined after trials. Such imperfections 
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introduce error in the system and reduce the accuracycof the system. 

Due to this, the probability of 'kill' is reduced considerably in practice. 

*4 Possibility of & Correction Technique 

Wien a digital computer is used alongwith the -weapon control 
system, the predicted future positions of the aircraft at regular time 
intervals can be stored to find the error in prediction at various time 
intervals. If the variation of -the error is uniform, these error values 
can be used to apply a correction to the predictor as explained belov. 

If the error in prediction of the range is DE^^ when the predicted future 
range is E^^ for various i, then DE^^ can be expressed as a function of 
E^^. llttien the coefficients of this function are kno-wn, for any predicted 
future range E^, the error in prediction can be computed. This is 

applied as a correction to E to reduce the error in prediction of . 

I3C -i-X 

This method of applying corrections is explained in the next two chapters. 

1.5 Outline of the thesis 

The following is done as part of this thesis works 

1 ) Simulation of the Weapon Control System on a 
Digital Computer. 

2) Development of an error correction technifine to improve 
the performance of the predictor and hence the overall 
system. 

5) System design of a special purpose Digital Computer to 
be used with -this Weapon Control System. 

The development of the correction technigue is arranged as part I 
and the system design of the Special Purpose Computer as part II. 



CHAPTER 2 


ERROR m EREDICTIOR 

As mentioned earlier, error in prediction might depend on many 
factors like the hearing, speed and eleTation of the aircraft and the wind 
conditions. This error is made up of two jarts , namely, 

1 ) Error in prediction of future position of the aircraft. 

2) Error in prediction of" the Shell trajectoiy and Shell 
position or the time of flight of the Shell. 

These are shown in Fig. 2.1 and will he explained helow. 

2.1 Total Error in Pcediction 

The Radar and the Gim are at the same point G, as assumed earlier. 

At the time, the aircraft is detected^ the CLOCK set to zero and the aircrafi;: 
position is indicated hy I. Depending on the present rates of the aircraft, I 
the predictor predicts that if a Shell is fired with an initial -velocity 
at an ele-vation of 1^, the Shell and the aircraft will meet at P after 
a time lapse of t^ (l ) where t^(l ) is the predicted time of flight of the 
Shell. The predicted future position, P and the time of flight of the 
Shell t^(l ) are stored. As the sampling rate is chosen to be 100 msec as j 

explained in a later section, the clock is incremented by 100 msec and the j 

future position and time of flight of the ^ell t^(2) are computed and | 

stored. This is repeated for every 100 msec. 

Ihen the clock eq.uals )» radar inputs are sampled. Let 
the aircraft be in position A, Hence PA gives the error in the prediction 





Flgi:i, 2«2 Shell Burst Pig. 2.1 s Pirrcr In Prediction 
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of the futixce position of the aircraft. At this time, the position of 
the Shell should he determined. But it is difficult to keep tracking 
a Shell with a Eadar. Also normally a Shell will not he fixed untill the 
correction function parameters are available . Hence the Shell trajectory 
and its position will he determined mathematically, from the aerodynamics 
and wind conditions. Let the position of the Shell, as determined by 
this mathematical model he S, Thus ES gives the error in the prediction 
of the Shell position. The total error in prediction is SA. 

When the Shell hursts at point S, the splinters from the Shell will 
he sprayed over an area X as shown in the Pig. 2.2. This area is called 
the destructive area of the Shell. If the body of the aircraft is anywhere 
inside this area, that is, if the areas X and Y in the figure intersect,, 
the aircraft will be hit. To achieve this objective, we will have to 
reduce the error in prediction to within + 5 yards. This means an accuracy 
of the order of 0,01^ at a range of 3OOOO yards. 

■We have assumed a sampling rate of 0.1 sec ^diich needs justification! 
Even an aircraft flying at a speed of 1 Mach will cover only about 55 yards 
in 0.1 sec time interval. Hence vhen an aircraft is in the range of | 

5000 yards to 3OOOO yards, samplings in steps of about 50 yards will he 
quite sufficient. Hence the sampling interval is chosen to he 0,1 sec. 

2.2 Simulation of the System 


This weapm control system is simulated on a Digital Computer 
and its performance studied. The predictor is represented by a set of 
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simultaneous equations. As mentioned earlier, a stable platform is 
assumed for the time being, hence the stabilizer is eliminated. How 
the problem is to simulate the aircraft (or the Badar) and the Shell 
movement. !Ihese are explained below. 

2.21 Simulation of Aircraft 

The initial position at which the aircraft is detected, is chosen 
at random with random number generators, to suit the typical cases that 
might occur in a real system. The speed and course, i.e. direction in 
which the aircraft moves, are also chosen. It is assmed that the aircraft 
maintains its initial speed and course. This assumption will be justified 
later. %enever there is a change in either the speed of the aircraft or 
its course, the system is initialised. 

If V is the speed of the aircraft, it moves a distance of v x dt 

during dt^ the sampling time (interval of 0,1 sec). Beferring to figures 

2.5 and 2.4, F gives the position of the aircraft after the time interval 

of 0.1 sec. H , S and T give the position of the aircraft at time t, 

P P ^ 

and B , S and T„ give the position of the aircraft at time (t+dt). 
and T^ are computed as shown below. Here, 

X»- Angle between the present bearing T^ and the course 
of the aircraft. This angle is called the angle of 
approach of the aircraft. 

ingle between the horizontal and the course of the aircraft. 
This is called the diving angle of the aircraft. 


and Y*- 
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The distance mcnred "by the aircraft in the time interval dt 

is greatly exaggerated' in figures 2»3 and 2.4. R„, S and T are 

i X s 

computed as shown below. 


Now as 0 is very small, 





II 

Sec Z' 

fji] 


= s 

p 

+ & 

[l.ij 

"f 

= T 

s 

+ Z 

[1.3J 

small. 




tan © 

11 




AP 




“* R _ AP 

P 



Vi:dt . 

eosX.Sin (Sp - y) 





= . OF = OB - ZB 

= Rp, Cos 0 “V. dt. CosX, Cos (S^ - Y) 


Also as Z f is very small, 

Z ' = tan Z ' 

GF 


OC 


* V. dt. SinX. CosY 




Deflection on horizontal plane, Z is given by, 
Z' 

^ ~ Cos S 

P 

V. dt. SlnX. CosY 
R^ . Cos 
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2.22 Simulation of the Shell Tra.iectory 

The Shell moves in super sonic speed and is influenced by 
the following forces* 

1) Acceleration due to gravity, g 

2) Shrasite drag due to air resistance 

3) Drag force due to shock waves and skin friction 

4) Wind effects 

, 1^ ) Gi^yX) i‘Co\>'0 Cy> IbJ ' 

Of these, the second and third are combined to get the total 
drag coefficient. Hie total drag force acts in a direction opposite to 
the direction of motion of the Shell. This can be represented as, 

S', tfC C .A. r5-43 

0. a 

where '^d*" Combined drag coefficient 

V *- Velocity of the Shell 

and A *- Area of cross-section of the Shell. 

The combined drag coefficient C^ is a function of the Shell velocity®' 

i.e., = f(v) 

can be experimentally determined for various' Shell velocities 
in a wind tunnel for any particular Shell. The variation of with V is 
given in Hg. 2.5* 

During the flight of the Shell, the velocity of the Shell vanies 
continuously. Hence would also vary with time and it will not remain 
constant. The velocity of the ^hell at any instant varies with the 
acceleration on the Shell. The acceleration (or retardation) depends on 
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the total drag force ishich in tijm depends on C^. But is a function 
of V, Thus a looping occurs as and V are inter-dependant. 

O-oil 


0*0 4 ^ 

a- ! I 

o-©3‘ ) 


/ ; 

l\ 

! \ 


0 - 01 *-^ 


0[ 


/ 


\ 


'Rr CmAchJ^* ^ 

- VotYtO. tton. Cd Tv 


To solve this problem tje Trill break the total time of fli^t 
of the Shell into a number of small equal intervals, as shorn in Fig- 2.6. 
During each small time interval, is assumed to remain constant. Having 
fixed C^, the total force on the Shell and hence the velocity of the Shell 
at the end of this time interval can be computed. For each time interval 
0^ will be different and is computed from the velocity of the Shell in the 
beginning of that time interval. 


Fow the number of such small intervals are to be chosen. To get 
this, for the same G-un elevation, time of fli^t and muzzle velocity, 
equal time intervals o£ 4,6,8,10,12 and 16 were assumed and the position 
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of the Shell at the end of the tine of flight was computed. It has "been 
found that there is little variation in the computed Shell position when 
the number of time integrals are more than 8. Hence the number of equal 
time intervals are chosen to be 8. 
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2.25 Error due to imnerfections in Simulation 

The mathematical model used to simulate the aircraft is an 
approximation to the real system and this introduces error in the solution. 
Similarly the analytical determination of the trajectory of the Shell also 
would deviate from the actual sitUfttlon. Hence the error computed 
as shown in section 2.1 between the Shell position S and the aircraft 




position. A (Refer to Pig.2.l) obtained from the simulation, program 
would comprise of the following* 

1 ) Error in the determination of the trajectory 
of the Shell. 

2) Error due to the non-idealness of the aircraft 
simulator. 

5) Error in prediction. 

It is felt that the first error can be reduced to a negligible 
value by collecting sufficient data regarding the Shell trajectories 
through a nnmber of experiments and trials under varrying circumstances. 

In a real vreapon control system, aircraft simulator will be absent and the 
aircraft positions are obtained from the Radar. Thus the error 2 ?hll be 
absent in a real system. This leaves the error in prediction as the major 
component of the error in the real system. 

Earlier it was stated that the predictor assumes that the aircraft 
maintains its speed and course. Thus the predicted position becomes inva.lid 
if the aircraft changes either the speed or its course. Tiiere are some 
predictors which assume a curved path that can be represented by a 
mathematical expression. Eut such predictors will fail if the aircraft 
takes a strai^t line course. Thus we have to compromise between these 
two different approaches to the solution of the fire control problem. In 
a real situation, the aircraft will pose more danger when it maintains 
a steady straight line coxrrse. Hence the predictor for the weapon control 
system under consideration is designed to tackle the aircraft effectively 
when the aircraft takes a steady strai^t line course. Thus the earlier 
assumption is justified. hmi ^ 
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COBRBCTIOH FIMCTION FOR QBE ERROR BT FREDICTIOR 


5 «’' Analysing the Error ig Fcediction 

In tile previous chapter, we have seen that the error in the 
system is mainly due to the error in prediction. To improve the 
performance of the system, the error in prediction is to be minimised. 

To develope some means of reducing this error, the variation of this 
error under various conditions are studied. The error might vary with 
the following main parameters* 

1) Present position of aircraft 

2) Speed of the aircraft 

3) Wind effects 

To take into account, all the possibilities that might occur in 
real situations, about fifty randomly selected typical cases have been 
analysed. These cases are tabulated in Appendix - I ■*. For each case 
the folloTfing variations of error in prediction with predicted future 
range is studied. 

1 ) DR Ys R^ ( R« Range ) 

2) DB Ys R^ (B i Bearing) 

5) DE Ys R^ ( E * Elevation) 

Some of these variations are shown in figures to 5»l(d). 

From these figures, the following points have been ob served * 

1) The variation of DB with R^ is mostly linear. 





1 ? 
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2 ) !Ilie variation of DR with. R^ is nonlinear. 

5 ) Ihe variation of DB with R^ is nonlinear. 

4 ) Hie variation of all the above mentioned errors with 
is smooth. 

These observations indicate the possibility of fitting curves 
to represent the variation of errors with R^ and extrapolate them to find 
the errors in future predictions. 


5*2 Curve Fitting 


The variation of DB with R^ is fairly linear as seen from the 


f3‘D 


jPi-g»3*l(a) to 5»'l(d). Thus a straight line of the form 

DB = A + , A . « E _ 
oil 

1 

can be fitted to represent the variation of DB with R^. For ten 
predicted future ranges E^.^ throu^ ^■flO' bearing DB^ through 

DB.^q are obtained. With these ten sets of data, the parameter A^ and 
are determined using linear least squares method, of curve fitting. Row, 
for any future predicted range R^^, probable error in bearing will be ■. 


DB as given by, 

X 


DB = A + A, . R 
X 0 1 fx 


£ 

Hence, if B is the predicted bearing, this is corrected by an amount DB 

X ..t 

Thus the an^ to which the Guns are to be trained before firing is given 
by the following equation* 


B = B - DB 

XC 2C X 

= B - A -A. ,E^ 

X o 1 fx 


r3-3j 


Bib variation, of DE with. and DR with are non— linBar* Thoy 
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can be conveniently represented by polynomials of the form, 


m = + C^. + C^.R^^ + 

m = R^ + D^.R^ + 


. . fS'^; 
. . Cb 'sJ 


To simplify the solution, only the first three terms of the 
polynomials are considered. Then the above equations will become. 


m . O^+C^.E^+C^.E^ 


DR . +1)24 fS'lJ 

These equations are linearizable as shown below * Rewriting the 


equations, 


where 


HE . 0„ + C^. E^ + C^. 

DR = D 4- D^ • R^ + D^ * R 

0 i I Z s 

\ - 4 


C^'SJ 


These equations for IE and DR are linear in R.* 3.nd R and their 

I s 


coefficients are computed as shomx helows 


A = DB 
o 


Ai - 


% R^ 

(Rfi - fif ) . db^ 


a, - (eg . DR - EA. IE) 

Q> « 

^ DP 

n (ER , EA - EE , AE) 

DP 


DE - C. * - C-. R 

1 I Z B 


( RE . DR - BA, AE) 
DP 


C^' "3 

fa* 1^3 

Cs- >3] 

(- 3 . I4J 
C 3- isj 



21 


(ER. Ril - ES . JiE) 

D. =; ■■■'■"_ 

1 DP 

C3-ik) 

®o = ^ - ^1 * ®f - ^2 \ 

[^•nj 

and DP = ( ER. DR - AE , AE) 


where DR = ^ (r^^; - R^). 

[3^qJ 

ER =;s: (E , - R ) . R . 

SI S SI 

[3-Zo3 

m = - r^ ) . de^ ■ 

1 S-2'J 

EE - 5 ),. BE. 

C h-iZj 


-0 


EA 

■= - 


) • 


-2hJ 

AE 

- - 


).R . 

SI 



dr” 

* Mean of 

DR 





1 Mean of 

DE 




\ 

1 Mean of R^ 




R 

s 

* Mean of 

R 

s 






A flow chart depictijag ''the correction technique is presented 


in fijj.3»2. Referring to the flow chart, when an aircraft is detected, the 
system is initialised and the real time CLOCK is set to zero, dhe Eadar 


inputs are sampled and the predictor computes the future position of the 


aircraft and the time of flight of the Shell. Ihe time of flight t^^ is 
added to the CLOCK and stored as t^(l ). 'Ehe trajectory of the Shell is 




Pig. 5,2- - PLOW CHART OF THE! mOi'rBm 
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determined and its position after t^^ seconds ^is computed as explained 
in section 2.2. This position in terms of range, tearing and elevation 
are stored as respectively. Now, after 100 msec, 

i.e when the CLOCK eq_uals 100 msec., the Inputs are sampled again and the ^ 
future position aad time of flight of the Shell, ■t^2 computed and stored as 
R^(2), 5^(2), and E^(2^. The CLOCK is added to t^g and stored as t^(2) 

This is repeated for every 100 msec time interval as we have chosen the 
sampling interval to be 100 msec. 

When the real time CLOCK equals t^(l), the radar inputs are 

sampled to find the actual position of the aircraft as given by E_, T 

P ° I 

and These are compared with R^(l), B^(l) and E^(l) respectively to [ 

get the errors LR (1), LB (l ) and LE (1). igain, the inputs are sampled 
when the CLOCK equals t|.(2) to get the errors DE(2), BB(2) and LE(2) and 
this process is repeated. 

A set of ten values of LR, LE and LB are collected for fitting 
curves between the following pairs of variables ? 

1) DR Vs R^ I 

2) LB Vs R^ 

5) DB Vs R^ ! 

These curves will be used as correction functions* Now onwards, for each 

input sampling the futxue positions R^Gc), B^(x) and Ej(x) as well as theirj 
corrected values and E^^(x) are computed and stored. If 

the aircraft is wilhin the firing range of the Oun, the firing can 
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commence with the corrected values of the future position i.e. (x), 

X o 

B^^(x) and 

As before, idienever CLOCK equals t^(x) the inputs are sampled 
to get the actual position of the aircraft at that time* These are given 
Rpj and S^. Now the errors in prediction of range, bearing and 
elevation with and without correction being applied are computed as 


indicated below. 

I)E(x) = E^(z) - 

= E^(x) - Sp 

DB(x) = B (x) - T r3- 

X o 

DE^Gc) = R^^(x) - E^ 

DE^(x) = 

BB^Cx) = B^^(x) - 


where, DE(x), nE(x) and BB(x) are error in prediction of range 
e elevation and bearing respectively without the correction 

being applied*. 

BE (x), BE (x) and BB (x) are the error in range, elevation 

C 0 c 

and bearing respectively with the application of correction 
technique* 

E„^(x), and B. (x) give ihe corrected future position. 

X C X c xc 

The Yariation of the error in prediction with and without the 
application of the correction techni<iue has been studied for fortjr eight 
typically chosen cases* These cases are tabulated in Appendix I# 
variation of the error in a few cases are already shown in Fig. 5 *1 (a ) to 
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The improvement in the performance of the predictor due to 


the application of the correction technique isc evaluated as indicated 


below* 


where 



...... I>F(x) 


DR (x) 
c 

^^(x) - ■ 

DB(x) 

[1 ■ 33i 

DB^(x) 

F^(x) = - 

...Mx)^ 

r 3 ' 34J 

BE (x) 
c 


(x)» P^(x) & improvement factors in range, bearing 

and elevation respectively. 


The variation of F^(x), and T’gCx) with R|.(x) for four 

typical cases are shown in Fig, 5.3(a) to 3»5(<i)* After studying these 
variations for the 48 oases given in Appendix I, the following conclusions 
are drawn s 

1 ) The error in the prediction of range is minimised 
(to a very small value ) after the application of the 
correction technique. It is ali^jays less than + 2 yards 
for slower aircrafts. For faster aircrafts DR^ is 
within + 5 yards in many cases and in a few cases 

DR is as high as + I5 yards. 

2 ) The improvement in the prediction of bearing 
elevation are much less compared to the improvement in 
the prediction of range (Refer Fig. 3»5(®‘) to 5»5(<l))« 
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>3) Soon after the correction coefficients are evaluated, 
the improvement in prediction is high tut as time 
passes on the improvement factor decreases appreciably. 

^. 4 ) spikes in the curves appear due to the variables 
DE or DB or BE cross the zero axis and has no 

G C C 

physical significance. 

5-4 further Improvement in Performance with Modifiei^ 

Ehe observations made under above indicate that the extrapol- 
ation of the correction function curve, with increasing time introduces 
larger disparity. To alleviate this and improve the performance of the I 
system a modified technique of curve fitting with latest set of data I 

points has been tried. 

This method of curve fitting is represented schematically * in | 

Fig, 3 . 4* This resembles a push down stack. Wienever N number of new | 

samples of error values are obtained, the ctirve fitting routine is entered 1 
and a smooth curve is fitted with the newest 10 samples. Thus the - i 

correction function coefficients are modified after every E number of | 

samples are obtained. This method has been tried with different F like 

I 

1, 10, 20, 40 and 80* It has been observed that there is ■ Tery little 
ixaproveinent in the performance of the predictor when the N is small 
compared to larger values of U* The results are shown in Elg« 5*5‘*(25‘)to(e)* 
A sample case with U as 100 has been chosen for investigation* It has 
been observed that the error in prediction with correction being applied, 
remains within the tolerable limits when the value of N is 100. Hence it 
is decided to have N as 100* 



o o 
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H new values enter 





Samples used for 
^ curve fitting 


Fig, 3 *4° Mo dified set cf data for curve fitting 



.■'-'R^(in thousand yards) 

Fig. 3.5 (a) Variation of F^ with N 


20 



30 




Pig. 3.5(c) 
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So far it has been assumed that the cxorve fitting routine is 
entered and the correction function coefficients are determined with ten 
samples of error values. This needs justification. The curve fitting 
has been tried with samples like 5, 10, I5, 20, 4O and 80 and the 
improvement factto in each case has been studied. The improvement factor 
for various number of samples is shown in ^ig. Fnom this, it is seen 

that having more than 10 samples has very, little effect on the improvement 
factor and it only adds to the complexity. Hence having only 10 samples 
for curve fitting is justifiable. 

5*5 A Critical Analysis 

In the sequal certain important aspects such as the cost of imple- ; 
mentation of this technique and its complexities, are considered. For 
these discussions, a special purpose computer for the Isapon Control System) 
is assumed to be - available. In Chapter 4 the description of such a Special) 
Purpose Computer will be given. ■ i 

3-»5’' Time Delay In Peteomiinina the Correction Function Coefficients 

The correction function coefficients are available only when the | 

I 

aircraft is at a range of about 20000 yards, though it was detected at a j 
range of 3OOOO yards. That is, when the aircraft is detected at 5COOO 
yards range, the present range is 3OOOO yards. For this, the predicted 
future range will be around 20000 yards when the aircraft speed is about 1 

'I 

600 knots. Hence only idien the aircraft comes within 20000 yards range, 
the first error value can be determined. As the curve fitting is done 
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when, ten such samples of error values are available, the aircraft will be 
within 20000 yards range when the correction function coefficients are 
determined. This much time delay may not be acceptaijlej To mininiize this 
delay the correction function coefficients are to be obtained much earlier : 
for which suitable methods should be provided. 

One method of getting over this difficulty is to assume a 
hypothetical Shell velocity idiich.is much higher than the act'ual Shell 
velocity. A Shell velocity of ten times the initial value has been assumed' 
in the beginning, to get the correction, function coefficients earlier. 

Now, if the present range is 3OOOO yards, the predicted future range would 
be around 2^000 yards, '^us the ten samples of error values required for 
curve fitting can be obtained \dien the aircraft is at a range of 29000 
yards itself and the correction function coefficients can be determined. | 
Once these coefficients are determined, future computations areldone with | 
the actual Shell velocity. 

Now we have to see #iether the use of these correction coefficients 
determined with the hypothetical Shell velocity, when applied to the 
computations with actual Shell velocity improve the performance of the 
system. After analysing the improvement factor with this method for about 1 
ten typical cases, it has been concluded that the improvement in the I 

performance of the predictor is only marginal. ?he improvement factor 
is not more than 2. Hence it is decided to t^e these correction function 1 
coefficients obtained in this method till the correction coefficients 
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can be determined with actual Shell ■velocity. For the case discussed just 
nov/, for the period, •fclae aircraft moves from 29000 yards to 20000 yards, the 
correction function coefficients determined -with the hypothetical Shell 
veloci'fcy will be used and when the aircraft is within 20000 3rards range, 
correction function coefficients determined with actual Shell veloci'ty will 
be used. 

5 •52 Cost of Implementation of •the Technique Vs Im'Frovement in Ferformance 

Ihen this correction technique is to be implemented in ’the existini 
system, this adds to the c'omplexi'ty of the system as given be lew J 

1 ) Additional storage space is required for the program. 

2) Additional da-fca storage for storing •the predicted 
future position and samples of error in prediction. 

This also increases the execution time of the programme. 

This simple technique is independant of the system and is 
applicable for any system using a predictor. Also even "when the aircraft 
is out of firing range, i,e during ■the initial tracking period, correction 
function coefficients are determined so that when the aircraft enters the 
firing range, it "will be possible to shoot do^wn the aircraft with better 

accuracy. ■ i 

As it has been seen earlier, the improvement factors 5^(x), 

F (x) and 'F:'’ (x:) are ‘.of the 1 order- .of i 10. '.Hence if a weapon control system 
b' e 

having jfo probability of ’kill' is chosen for the implementation of the 
technique developed, the probability of 'kill’ of the modified system will I 
be above 40^. Tn general, 
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p =i-(i - p r h' 

where P^s probability of 'kill' of the modified system 
PqJ Probability of 'kill' of the original system 
and n * PnproTement factor of the correction technique 

^uch an improvement in the system performance Is of great 
advantage in Weapon Control Systems. Even an improvement factor, of 2 will 
justify the cost of implementing the correction technique in Weapon Control 
Sys terns • 
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A SEBCHL PUBBOSE COMPOTBR. FOR lEiPOH CONTROL 

In this chapter a special purpose Digital Computer (used with 
a weapon control system) will he discussed. In military applications j 
there are stringent conditions on the factors like power req^uirements, 
space requirements, reliability and temperatixre range in which the system 
has to operate in addition to techniques for the removal of the dissipated 
heat. 

Digital systems using vacuum tubes had a number of problems like 
heat dissipation>size and reliability. The mean time between failures 
(MTBP) of such systems were only few hours. Developments in the semiconductoi 
devices and the Integrated circuits (ic) technology have made it feasible to ; 
desifjn smaller and less costlier but more reliable systems. The high 
packing density of iCs, has further reduced the power consumption of such 
special purpose digital systems. Da view of this, the system described 
below will have Integrated Circuits design. 

4.1 Dse of the ICs in Digital Systems | 

Ihen small units like gates and units such as flip flops which use 
less than 9 gates are formed on a single monolithic chip, it is called | 

Small Scale Integration (SSl)^ Ihereas in Ifedium Scale Integration (tISl), 
logic sub assemblies such as counters, shift registers, decoders, memory 
elements, full adders and look ahead carry generators vhloh. use less than 10( 
gstos srG foriuod on a singl© monolithic chip* More complex circixits which : 
need more than 100 gates such ^Anslog to Digital (iy^) convBrter and 
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Digital to iinalog (d/a) converters are made in a single monoli'thic chip 
in large Scale Integration (LSI). lotting larger sub-systems on a single 
chip has the following advantages? 

1 ) Reduced Space and weight 

2) Higher reliability due to reduction in external inter- 
connections. 

3) Easy maintenance due to the minimization of test points. 

4) lass cost per Unit. 

Designing systems with IC sub-assemblies is much simpler. Further 
there is reduction in inventory and spare parts stocking. Ihe failure rate 
of individual gates are of the order of 3x10 per hour. Hence the tflBF 
for a computer using IC sub assemblies will be of the order of few thousand 
hotors. Hence the computer described here will use ICs . 

In Chapter 3j S' clocking rate of 100 msec has been considered to be 
sufficient for the weapon control system under discussion. Hovjever a 
provision for flexibility in the system is required to deal with the improved 
speed in future aircrafts and to perform additional functions as they mas'" 
be required. Hence it is better to have a speed improvement factor of 10. 
Tlius the special purpose computer under discussion should be able to solve 
the fire control problem within an optimum time of 10 msec. Due to the 
practical considerations, the system should have minimum power dissipation 
and the cost should be reasonable. 

Integrated Circuits use many different logic types. Some of them 


are given below* 
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1 ) Direct coupled transistor logic - DCTL 

2) Resistor - transistor logic - RTL 

5) Resistor - Capacitor-transistor logic - RGTL 

4) Diode - transistor logic - DTL 

5) Transistor-transistor logic - TTL 

6 ) Emitter coupled logic - ECL 

Of these, ECL is the fastest circuit, having switching time of less 
than 5 j&sec with power consumption of atout 50 mw per gate. I’/hereas IfL 
family has a wide variety of gates for different applications. The standard 
sorioa TTL gates have typical propagation delay of 9 to I5 nsec and poY;or 
disoiiyition of about 1Qmw per gate. High speed TTL series have switching 
speed of about 5 usee but the power dissipation is about 55 ^iw per gate. 

Low powjr TTL aeries gates have power dissipation as low as 1 rmt per gate 
but the propsigation delay is about 35 nsec. Of all the 10 logic types, TTL 
gives the highest performance to cost ratio and among the various TTL 
series, the low power TTL series has the best speed-power product. But the 
low power TTL series are many times costlier than standard series TTL. Hence 
the low powr TTL series units or sub assemblies can be used in systems 
whore the cost factor is not very important compared to power requirements. 

In the system under consideration, standard series TTL circuits will be used 
as the increase in power requirements due to this, is reasonable, -un the 
computor has to operate in the temperature range of -55°G to +125°C, logic 
sub assemblies designed for the operation in this temperature range, such 
as the 5400 series (of the Texas Instruments Ihc.), will be used. The 
flow diagram of the proposed system is given in Fig.4*1* , 














40 


4*2 of th© instruction Set ' 

Considering the various calculations involved in the solution of 
the fire control problem, a basic set of instructions given below are 
proposed in the first round. 

1 ) Loading a register from memory 

2) Addition 

3) Subtraction 

4) Multiplication 

5) Division 

6) Storing the contents of a register in memory 

7) Test for positivity 
s) Complementation 

9) Sampling the input 
10) Outputting (the computed results) 

Fi^rthei" provisions are to be made for generating functions, such as. 

1 ) SINE (X) 

2) cesiNB (x) 

3) AECBIBE (X) 

4) SQUARE ROOT (x) 

Biis can be done in two mjs. The first is a software approach 
whore subroutines will be stored to realise these functions. The second 
approach is to store these functions in microprogrammed memory and realise 
them hardwarewis(^. First, the software approach will be tried to see whether 
this will suit the requirements of the system under discussion. Having 
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microprogrammed memory adds to the complexity of the system. Hence the 
hard..^.re approach .rill be considered only if the software approach is not 
sxiitQ<bl© for tho systGni tmdGr considoratioti* 

Tiie following mnemonics are chosen for further discussions. 


Maemonics s 


Operation 


1 ) 

2) J^S 


Load the Accumulator Eegister with the contents of the 
memory location addressed. 

Add the contents of the addressed memory location to the 
contents of the Accumulator (Aoo ) and leave the result 


5 ) SB I 

4 ) Mils 

5 ) m. 

6 ) ST£ 

7) Sh 

a) cs* 
9) TP*, 


in the Acc. 

Subtract the contents of the addressed memory location from 
the contents of the Acc and leave the result in the Acc. 
Multiply the contents of the Acc by the contents of the 
addressed memory location and leave the result in the Acc. 
Divide the contents of the Acc by the contents of the 
addressed memory location and leave the result in the dec. 
Store the contents of the Acc in the memory location addressed 
Ebcecute the addressed subroutine with the contents of the 
Acc as the argment. 

Complement the contents of the Acc. 

If contents of the accumulator is positive transfer to 
the location whose displacement is given in bits 4 through 
11 of la (Instruction Eegister), else contiirue. 
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10) 

BN* 

Hetum to the routine from where a subroutine 

was called. 


11 ) 

SPs 

Sample the input indicated by bits 8 through 1 1 of 

and store it in the Acc. 

IB 

12) 

OUs 

Output thG couteuts of the Acc, through D/A converter, to 

the output terminal indicated by bits 9 through 11 of IH. 


4*5 Mamory Bequirements 

With the instruction set given above, the fire control problem is 
programmed and is given in Appendix-II. Stom this, the following can be 
observed* 

1) ybraoxy required for main program = 276 locations 

2) Memory required for SIl/cOS Subroutine = 45 locations 
5) Memory required for SQET Subroutine = 21 locations 

4) Ifemory required for AESIN Stibroutine = 37 locations 

Thus total memory requirements of the program = 379 locations 

This special purpose computer operates essentially on a fixed 
program. From Appendix II, the memory requirements for data storage is 
91. Hence 5I2 words of memory will be sufficient for the computer. 

4*51 Tata Format* 

The ratio of the largest number to the smallest number to be 

15 

handled by the computer under discussion is about 10 . If fixed point 

arithmetic is to be used, about 44 bits word would be required to store the 
data. Hence 48 bits word length memory would be sufficient. 
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liereas if the 48 bit word length is considered xmeconomical, 
floating point arithnatic can be used. In floating point arithmetic a 
v".lx; bit signed exponent can handle numbers in the range of 2— which is 
sufficient for the system under consideration. Further, the allowable error 
in the system is + 3 yards (see section 2.1 ). Hiis is 0.01 percent of 
30000 yards which is the maximum range of interest. A choice for the word 
length as 24 bits with six bits for exponent and 18 bits for signed mantissa 
v/ould satisfy such requirements. Thus the use of floating point arithmstic 
would roquiro ninety one 24-bit word length memoiy locations. This saves 
ninety one 24-bit word length memory locations. But the use of floating point 
arithiifs tic adds to the complexity of the system. Further, provision for 
convtirting floating point number into fixed point representation and vice 
versa are to be provided. Considering these, fixed point arithmetic will be 
ohoson for the oominiter, allowing for tolerance of minor error in 

calculation. 

4.32 Instruction format 

Tho total memory requirement is less than '^00 words. Hence 9 
bits aro required for the address part of the instruction. For the t^^relve 
oferation codes (opcodes), four bits will be required and the instruction word 
length becomes I3 bits. A different organization is considered to use 12 
bit words fur storing the instructions. This will allow only 3 bits for 
opcode. With throe bits, eight combinations are possible. As some of the 
opcodes do not need an address part, the eighth combination, namely, 111 
is used to indicate such opcodes. The opcode table is shown below: 
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1 ) ID T 00 

2 ) iJ) 000 

3 ) SB 001 

4) MU 010 

5) DV oil 

6) ST 101 

7 ) EN 110 

8) TP 1111 

9) SP 111000 

10) OU 111001 

11) GS 111010 

1?.) EN 111011 

Iki it is seen from the program given in Appendix II, the 
diaplncomont for tho conditional transfer instruction is al-ways less than 
huri'ln.'d locotione* Hence tho signed displacement for this instruction will 
need 8 hits namely bits 4 through 11 of I.R. Hence the four bit opcode 
1111 is chosen for this instruction. Hence the first foxur bits of the 
^'U!i.'unin/' four opcodes, namely, 3P, OU, CS and EU should be 1110. Hence 
tfw-'ir I'jxu'dfjB are chosen to have 6 bits as indicated in the opcode table. 
Ibrthi.'r, inntruotion SP will use bits 8 through 11 to indicate the input 
tonninal to be chosen and OU instruction will use bits 9 throtigh 11 to 
indlcnte the output terminal to be chosen. 

4*4 Inout/Output Units 

Tlio computer under discussion will be executing a fixed program. 
Qnoo tho’ progrfwi is entered into the computer, it requires no modification. 
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In order to mijaiinize the complexity of the computer it is decided not to 
have complex Ihput/Output devices. However an operator's console is 
conaidcred osaontial. Ihe program can be stored in magnetic cards. The 
data can conveniently be stored in IC memory. 

Ihe computer will have Analog/iigital and Digital/hnalog converters. 
Ihese will bo used in multiplexer mode as there are thirtem inputs to be 
sampled and six computed results are to be given as output as indicated below! 


IniJut 1 

M 

Input 2 

« 

E 

Input 3 

Si 

input 4 

\ 

Input 5 

R 

P 

Input 6 

S 

P 

Input 7 

T 

8 

Input 8 

S 

Inifut 9 

DTD 

Ihinit 10 

E 

sp 

inimt 11 

B 

wr 

Input 12 

W 

Iniiut 15 

E ' 

Output 1 

T 

S 

Output 2 

E 

g 

Output 3 


Output 4 

R 

gz 

OutiJUt 5 

®ld 

Output 6 
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These samplings are done in time interval of 100 nsec as 
mentioned in chapter n, Itovision will be made for raising an alarm if 
cither of the two converter outputs remain constant for an amount of 
preset time, say, for more than 5 msec. This would indicate a possible 
failure of the system. 

4.5 Arithmetic Ihiits 

The following arithmetic units are required for this computers 
1 ) One 48-bit fixed point full adder 

2) One 48-bit fixed point multiplier 
5) One 48-bit fixed point divider. 

On fi s'Lngle IC chip, fixed point adders are available. These 

oircxiitf! nix.‘ oaimble of performing parallel addition and subtraction. Using 

ru'ixriited addition or subtraction combined with shifting right or left, 

5 

luultipliontion and division algorithms can be realized . 

The following set of registers will be required in the Central 

fniousoing Unit (CPU), 

1) T\vo 48-bit Accumulator Registers, Acc 

2) Qno 46 bit Multiplier Quotient Register, 

3) One 12 bit Instruction Register, IR. 

4) One 9 bit Location Counter, LC. 

5) Qno 9 bit Register to store the contents of LG 
while calling a subroutine, 'G* Register. 

6) One 48 bit Register to hold the addend during 
addition and subtraction, ’B' Register. 
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?) One 9“'bit Memory Mdress Register, MR 
B) Qao 48-bit Memory Buffer Register, MBR 
9) Quo 40-bit Input/Output Buffer Register, l/O Buffer 
*1' • ^ Ci-ntrcl Circuits 

Tliero are two basic cycles in the computer operation. They ares 
1 ) Inatrtjotion fetch cycle or I cycle 
2) Exocutlon cycle or E cycle 

During the I-cyole, an instruction is read from the memory and 
lu’^nicilit t'l tho IR* fflhiis Involves the following sequence of operations. 

'I’1 1 Shift the contents of the LC to MdR 
T?i Give a read pulse and iivait till the reading 
is complete. 

1’5* Shift the contents of MBR to IR. 

The circuit shown in Fig, 4.2a will generate the timing sequence 
pulnr.'s for instruction fetch. The timing pulses are shown in Pig.4.2b. 
in instruction fetch is initiated under the following conditions s 

1) Either the start signal is ’ON’ or the instruction 
complete signal occurs. 

2) Stop or Halt signal is ’OFF’. 

5) The above two conditions coincide with a cloclc pjulse. 

The following assumptions are made and justified: 

1) ’n' is the number of basic units of time or clock time , 

the I cycle takes. 
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2) The time between two clock pulses is greater than 
the rise time of the flip flop, FFI . This is true 
since the clock time is I50 nsec as chosen in a later 
section whereas the rise time of the flip flop-is around 
20 nsec. 

5) Bie clock pulse duration is less than the flip flop 
rint; time. Hence the clock pulse duration is chosen to 
be less than 20 nsec. 

4) n.ip flop rise time is greater than the gate delay-' time. 
'.Hiis is true as the gate delay is about 10 nsec. 

5) Start pulse coincides with a clock pulse. liiie is true 

as the clock time is I5O nsec and the start pulse duration 
would bo more than this. 

On completion of the I-cyole, the E-cycle will commence. During 
E-oycli.', first tlio opcode is decoded with the circuit shown in Fig.4.5. The 
c'luonce of micro operations to realize various opcodes are as shown beloy/J 

1 ) ID~100« Doad accumulator with the contents of addressed 

memory location. 

1) Shift the contents of the address part of IE to MIR. 
ii) Give the read pulse and wait till the reading is over, 
lii) Shift the contents of MER to Acc. 
iv) Switch the 'operation complete' signal 'OH'. 

2) iiD-OOOi Add the contents of the addressed memory location 


to the Acc. 
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i) Shift contents of address part of IR to mar . 

ii) Give a read pulse and Tisait till the reading is complete, 

iii) Shift the contents of MBE to 'B’ register. 

iv) Connect the Acc and 'B' register to the hardiware 
adder and start i^ddition Cycle, 

v) ii-t the end of addition cycle switch the operation 
complete signal 'OUT', 

5 ) S 1 Subtract the contents of the addressed memory 
location from the contents of the Aoc. 

i) Shift address part of IR to MA.R. 

ii) Give a read pulse and wait till the reading is complete, 

iii) Shift the contents of MBR to 'B* register, 

iv) 1’® complement the contents of 'B' register, 

v) Connect the 'B' and Acc registers to the adder and 
start addition cycle. 

vi) Tlien the addition is complete, switch the operation 
complete signal "CN". 

4 ) MU-010 1 Multiply the contents of the Acc by the contents of 
the addressed memory locations 

i) Shift address part of IR to MAR. 

ii) Give a read pulse and wait till the reading is complete. 

iii) Shift the contents of MBR to 'B' register. 

iv ) Connect the ' B ’ and Acc registers to the Multiplier and 

start multiplication* ; 

v) Ihen the multiplication is complete, switch the operation 
complete signal 


I f ! T k’AklPim I 
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5) P . lrQl. , 3 ,,^ Divide the contents of the Acc by the contents of 
the addressed memory location, 

i) Hhift address part of IR to MiR. 

ii) Clive a read pulse and wait till the reading is complete, 
xii) Shift the contents of MBR to 'B* register, 
iv) Connect the 'B’ and Acc registers to the Divider and 
start division. 

v) 'Miun the division is complete, switch the operation 
complete signal 'OR*. 

C'1'~1()1 1 Store the contents of the Acc. 

i) Shift the address part of IR to MiR, 

ii) Shift contents of Acc to IBR. 

iii) Give a write pulse and wait till the writing is complete, 

iv) Switch the operation complete signal 'OR'. 

7) S!K~110« Execute the called subroutine with the contents 
of the Acc as argument, 
i) Shift the LG contents to the 'C' register, 
ii) Shift the address part of IR to LC. 

iii) Switch the operation complete signal 'OR'. 

Q) S£dill* convention! Sign hit 0 — positive 

1 — negative 

Chock the sign bit of Acc,. If it is ”1 give end of 
operation signal. 

It the sign bit is Or 

i) Add the contents of bits 4 through 1 1 of IR to ID 
ii) Switch the end of operation signal 'OR'. 
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^ 1,9,9,Q.!.. Sample the input referred 

i) Boluct the input terminal number as given by bits 8 through 
11 of m and connect it to the A/d converter, llien the 
i/D conversion is over, iiifi output of A/D converter 

iri the Acc. 

li) Switch the end of operation signal 
Oilf 111001 1 Output the result referred. 

i) Ofinncct the Acc to the D/A converter. Connect the output of 
d/a otsnverter to the output terminal indicated by bits 9 
thfou-'h 11 of IR. 

ii) TOuin tile conyersion is over switch the operation complete 
oignul 'ON*. 

11 ) GS> 1110101 Change the sign of the contents of Acc. 

i) 1*® complement the contents of Acc. 
li) Switch the operation complete signal 'CUT*. 

12 ) BI 1 1 1 01 1 1 Betum to the main routine 

i) Shift contents of 'C* register to LC. 
li) Switch the operation complete signal *0N*. 

■Thus, a sub set of the following set of micro-operations 

have been uoed in all the micro programs given above. 

M1* Shift contents of ID to MAE 

B/Ei Shift address part of IB to MAE. 

151 Give read pulse. 

M 41 Shift contents of MBE to Acc. 

M 51 Shift contents of MBE to *B* register. 
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M 6 i Shift contents of MBR to IR. 

M7» 1 complcniont Acc, 

M 8 i 1 '® complement register. 

M9» Mdition cycle. 

MIOi Multiplication cycle. 

Ml 1 » Bivision cycle * 

M12* Shift oontonts of Acc to MBR. 

M 15 » Give a write inilse 

MI 41 Shift contents of IC to *G' register. 

MI 51 Shift address part of IR to LG. 

M16* Add bits 4 throng 11 of IR to contents of LC. 

MI 7 * Shift contents of *0’ register to liC. 

In addition to these micro- ope rations, the following 
Mnoro-o|)u rations have been nsedt 
MAC1 I k/h conversion. 

MAC2I B/A conversion. 

A gate level design of circnits for these Micro and 
Mii.cro-up..Tritiono are given in Appendix- HI. Ikecution of each opcode 
involves selection of a set of these micro and macro-operations. This is 
indicated in table 4«1 • ^ 0 ^ selecting these set of micro and macro operations 

eithur micro programmed memory or conventional circuits can be used. 
Micro/i&oro sequence generators (MSG ) use conventional circuits- Ihen an 
opcode is decoded, the corresponding MSG gets a pulse and this MSG 
generates the necessary pulses to choose the micro and macro operaf n 
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to Ik; cxncut.'cl. 4 level design of the liBGs is given in Appendix I?. 
Tilt) tijiiin^; pnlneo fjonerated by MSGJ to execute the opcode DY, is given 
boliW' net nn itxfMiiplc:: 


1 

ChOCK II 






It is observed from the list of Micro-operations, that many of 


thu inioro-operations involve shifting of contents of one register to 
antithor rogieter. Ihe time required for such operations will be chosen 
no tho baoic tijne of the computer. Biis operation will involve typically 
two IuvoIb of gate delay and the time for a J-K flip flop to toggle. 

TOum J-I flip flops designed for a toggle rate of 20MHz are used as with 
stantord series TIL gates, the total tire delay would he about 100 nsec, 
i olook mto of 150 hseo will satisfy the reluirement. 
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4*^^ for .executlo A of ,_the program and power 

'I'he number of times each instruction is executed and the time 
for E-cycle of various instructions, are tabulated in table 4.2. iinalog/ 

Digital and DigitalAnalog converters with 5OO KHz parallel word rate are 

5 

osnuiiiod . The optimum time to solve the fire control problem is 10 msec 
ao (t.'oicli.'tl in Section 4*1* '^ftthin this IQmsec, the following should be doneJ 

1 ) Sample all 1 3 inputs with A/D converters . 

2) Compute the output variables from the input variables. 

3) Convert the output into Analog signals with D/il converters. 

Fi'Djii table 4*2, it is seen that the execution time excluding 
I-ojolt^is aiJ])roxiinately 4 msec. 5he total number of instruction fetch in 
one iteration Is about 1600. Ihus one I-cycle can take a maximum time of 
a’bi'Ui; 3*7ytiSQc. Conventional core memories with 2 yasec cycle time would 
very liiuoh ocrve the purpose. 

4*9 AdvantUitjos of the proposed computer ». 

The special purpose computer described here has the following 
.'ulvantagoo when compared with one of the existing Electromechanical computers 
for uidvlng the fire control problem. 

1) Ifaintenance < The system described here, has logical 

sub-assemblies as basic unit. Hence maintenance becomes 
comparatively easy. Also this effects a reduction in inventory 
and spare parts stocking. 

2) Hfi liability t As explained in section 4. 1, the mean time between 

failure for the proposed system is of the order of few^ ^ 
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tliousand hOTjrs and hence the system is much more reliable 
than the existing system. 

5 ) In the existing system, there is lag in each servo 

unit. Due to this, the entire system becomes sluggish, 
l^ioreas in the proposed new system, the time for computing 
the results is less than 10 msec. Even this delay can be 
taken care off by applying a correction term to iiie predicted 
time of flight of the Shell, t^. 

4 ) Space * Iho electromechanical system will require a big room 
to install all the units. But the proposed new system can be 
kept in a small compact xmit (apart from, its 512 word memory 
and operator's console). 

5 ) Power consumption* In the electro-mechanical system using 
vacuum tubes, power dissipation and the removal of disspated 


heat pose a hig problem. Ihereas the proposed new system 
consumes less than 3 ^ watts power as can be seen from table 
4 .5.Ihia)L* less than 5 percent of the power requirements of 
the existing system. Further ICs can operate satisfactorily 
in the temperature range of -55°*^ "te +125 ^ lAdiich is good 
enough for the weapon control system. Cooling the new 
system is very simple as the amount of heat to be removed 


is very small. 

6 ) npATAtioual urooedurel Once the new system is ins tailed , it 

needs very little modifications. Testing is comparatively 

simple. The new system does not need constant attention of 
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8 bit, 52 word IC memory 1 9 


?, in])Ut, AND Gab 


5 Ini'Xit And iiiate 
4 Input ADD foito 
b ,to()Ut AMD 


Power 
per imil 

Total 

! Power 

1 

; vnW 

240 

1 2160 

10 

! 6680 
i 



9 input Oii 1 

1 I'Ui 16 injnxt aclocttv 1 






Lt ' k ,'Uiund carry genorat orj 4 
4 bit binary counters 11 


12 bit A/D Converter 


10 bit d/a converter 


Shift Registers- 16 bit j 15 


300 


500 


200 


TOTAL 


6.98 


5.85 


LlOQ. 

50.0 
3. 81 


35.0 


21.6 


52.88 


20.84 


33.66 


248.0 


282.5 


52.56 


1404.47 
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ttu; oporutor* Biere is no need either to tune the systesi 

‘vr to halojico the amplifier output as req^uired in the existing 
li.yiiti.'m, 

7 ) cost of the proposed system is expected to be arotind 

Rs.20,000; which is comparable to the cost of the existing system, 
due to the added advantages of the new system. 

^0 feotinnical mrtst There is no mechanical moving parts in the 
jin ij, list;, l syMtein. Yvhereas the existing system has wear and tear 
i!ui: ti) nu.'chaiiioal moving parts. ■ 
y) OthciT 0 1 ns idorntiona t In the proposed system with extended 
lusiraory, the error correftion technique as developed, in chapter 
3 can ho implemented to improve the performance of the weapon 
control system. 



CffilPTER ^ 

■^^^NCLUSIONS MD SCOEB fob FURTOTiIR 170RK 

X. prodiotoj? co3rr6ction techniq.ue has been developed and found 
to wonk StLtiBfactonily fon fonty eight typically chosen cases, which 
include most of the possibilities that occur in real situations. Ihey 
rt;juH'!!unt ni.fflhinntione of the following situationsi 

XI I Xlrcraft detected at long range. 

iX’i Aircraft detected at close range. 

hi I iippruaching Aircraft. 

B2 f Aircrvaf t moving away. 

Oil Ehster aircraft. 

02 f Slow aircraft. 

D1 I Low flying aircraft. 

DP* Aircraft at high altitude and diving down. 

Aircraft at high altitude and not diving down. 

Ell Calm atmosphere. 

E2I Moderately high wind speed. 

Tiio sjjfjoial purpose computer described here has many advantage 
over tho oxintlng oystam like, better speed, easy maintenance, better 
roliability, simple oix^rational procedure , lesser power cousumption and 

requires lessor space (see section 4.9). ^ logical organization for this 

ojjfcial purpose computer has been discussed. There is scope fo 
work in tfe® following lines! 
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1 ) tiufficient data regarding the course an aircraft takes 

under various real time situations can he collected and stored, 
'ftiia can ho used to account for the pilot's actions. 

2) 'Pho logical design of the proposed computer can be tested by 
simulation on an existing Digital Computer. 

5) Hie proposed computer can be fabricated and tested under 
the actual field conditions with a weapon control system. 

4) proposed computer can be modified to have extended memory 
find the predictor correction technic[ue as developed in Chapter 3 
can bo Implenionted. 

*■}) it. mmlifiod version of the proposed computer with extended 
facilities can be situated in a centralised place and used to 


control many weapon systems 


APEEUDK I 

hm . . OF TYPICAL CASES AWALYRTHT) 

'i’ll., fnlluwing is the list of 48 typically chosen cases which 
nnl. 


(Y-.rd 




(ibg . ) 

(Dog 

*" i 

(Deg.) 


(Kents, 

1 (KhDtsP (Dog. ) 


6076. y 

5 ^ 000.0 

50000.^ 


300.0 8.5 195.7 


43,6 180,0 0.0 300.0 24.6 104.3 


175,8 180,0 0.0 600.0 18.8 59.3 


600.0 38.4 257*0 


54.6 


280.5 


0.0 0.0 ! 300.0 

I 


C.oj 0,0 300.0 29.2 


,29.0 

280.1 

19,6 

158.1 


66fi6.4 


72.8 


296.8 


48.9 


351.8 


500.6 


15*5: 


45.8 


180.0 0.0 


0,0 


0,0 0.0 


180.0 -5*9 


600.0 


600.0 20.8 , 355.1 


300.0 9.2 320.6 


155.0 


12.4 552 . 4 


600.0 29.9 280.1 


15.8 21.1 


177*7 


Contd# • * 
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N'i . 

i. i 

li 

— 

T 

n 

-IsiiiL 

X 

.(De«.) 

Y 

(Deg.) 

Y 

(Saots ) 

w 

(fiacts) 

B' 

wr 

(De^ « ) 

19 

»•' -Hfc . . •* - IT 

8.0 

25.6 

180.0 

-4.9 

600.0 

13.6 

36.5 


•.117,'.* 

•5U( ;Ot),n 

fi.O 

22.7 

180.0 

-5.6 

600,0 

51.5 

153 -7 


H.O 

540.9 

0.0 

-5.6 

500.0 

14.8 

324.4 

?? 

50000.1) 

n.o 

160.0 

0.0 

♦4 

300.0 

35.0 

31 8.9 

2 ^ 

^0( !00,0 

a.o 

37.0 

0.0 

"5.6 

600.0 

16.3 

180.0 

2 .\ 

500f)r),n 

tun 

286.8 

0.0 

-5.5 

600.0 

35.0 

51 »6 

23 

6,079 .‘! 

2 »o 

145.6 

180,0 

0.0 

300.0 

13.0 

133.6 

i>r> 

1 

5 f,r{ .7 

2,0 

3.5 

180.0 

0.0 

300.0 

30.0 

206.1 

27 [ 

4055,5 : 

2.0 

158.0 

180.0 

0.0 

600.0 

17.9 

45-2 



............. 

5104.5 

2.0 

202.2 

180.0 

I 0.0 

600.0 

37.2 

55.0 

i 

?9 1 

17105,5 

2.0 

127.1 

0.0 

o.p 

300.0 

9.3 

I 186. 9 : 

50 

19204.1 1 

2,0 

21.3 

0,0 

0.0 

300.0 

I 25.8 

' 248.8 

31 

10701,0 



2.0 1 

155.3 

0.0 

0,0 

600.0 

: 13.4 ; 

1 10.7 

52 

15195.0 

2.0 

225.1 

0.0 

0.0 

600.0 

: 31*2 

1 

:-534.5 ; 
L -_4 


Contd , . . . 
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pRaiRAfvi 

TO SOiVE THE FIRE CONTROL PROBIHV! 


Tn lUilVi: 

tht! fii 

•0 control problem with the special purpose 

oomimtor 

di'-acribed 

in Chapter 4» the following 

: program is usedJ 




mm moa^M 


SP 

IHPirp 1 


SP 

INPUT 10 

S'P 

DT 


ST 

B 

SP . 

SP 

IIlillT R 


SP 

INPUT 11 

S'P 

h 


ST 

B 

wr 

SP 

INPUT 5 


SP 

INPUT 12 

Sti' 

§ 


ST 

17 

SP 

ttiPU'J' 4 


SP 

INPUT 15 

S'P 



ST 

E' 

q. 

SP 

W PUT 5 


ID 

R 

Sf 

R 

J' 

i:!piiT c 


MU 

dt 

SP 


AD 

R 

S'!' 

s 


ST 

R 


J* 




SP 

liiPllT 7 


MU 

at 

ST 

T 

t) 


AD 

R 

P 

SP 

urniT fj 


ST 

R 

P 

ST 

s 


ID 


SP 

IMFII'P 9 


MU 

at 

ST 

pro 


AD 

s 


Contd,... 

(CantittWjd in next column) 
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Rf 

i* 

hi 

m 

Sin S 

LlJ 

\ 

ST 

Acc2 

r.ii 1 

tit 

ID 

# 

S 

iiD 

N 

MU 

Cos S 

ST 

''i 

SB 

Acc2 

m 

a 

!• 

st' 

S, 

d 

m:! 

sra 

ID 

§ 

ST 

Dill s 

1' 

MU 

S 



ST 

Ago 2 


D 

I* 




ID 

1 

mi 

COS 


1 

ST 

Ci-o S 

MU 



i> 



LD 

s 

AD 

Ago 2 

BN 

s.m 

EM 

SQRT 

ST 

Sin S 

ST 

w 

i.W 

S 

ID 

^1 

affi 


MU 

Sia S 

a i 

lAiC 1;^!^ 


J 

LT) 

s 

m 

Cos S, 

i'ili 

cos 

ST 

C 

ST 

cos s 

MU 

h 

MU 

K 

ST 

¥ 


1 


* 


ii.cc2 

I® 

E 

iJ) 



» , 

MU 

CT1 

ST 

®ld 





MU 

w 

m 

K 




(Ctmtinuod in. the next ooltmin) 


Gontd 
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m 

E 

F 

MU 

CT2 

ST 

# 

w 

m 

CT5 

MU 

B ln C 

ST 

wc 

ST 

iiOc2 

lA LD 

R, 

LD 

§ 

ST 

"fo 

MU 

i 

C 

LD 

R 

i.D 

iHU’,',* 

lU 


ST 


AD 

R 

P 

LI) 

* 

« 

ST 

Acc2 

MU 

GtniC 

MU 

Acc2 

SB 


ST 

Acc2 

ST 

# • 

S 

ID 

E w 
P 

Ld) 


MU 



E^Sin K 



VI 


X 

ST 

E w 

MU 

R Sin D 

I 

F 

AD 

Acc2 

wUl 

w 

m 

SQRT 

ST 

li 


y) 

ei 

r 

ST 

R, 


ID 

RjSin D 

wu 

cos 

DV 


ST 

Cos Sji 

ST 

Sin D 

LD 




f 

M 

ARSM 

Mii 

sw 

ST 

D 

ST 

Sin 

MU 

Cos C 

LD 

m 


a 


ad 

(Continued in the next column^ 


Gontd* • • 



70 


ST 


ST 

Acc2 

IJ) 

U'll'l 

MU 

CT7 

Mil 

‘S- 

EN 

AES 11 

AD 

CTU 

MU 

CT8 

MU 

m 

ST 

TE 

ST 

Act;? 

m 

sm 

liD 

IHvr* 

MU 

CT9 

S:i^ 

r[) 

At 

ST 


S'lS 

DTD 

LD 


ST 

|%^» 

SB 

"fo 

EN 

SID 

TP 

PS 

ST 

Sin Uwf' 

CS 

« » • 

ID 

i'WK 

PS SB 

CT10 

E« 

(iOS 

TP 

u 

GT 

iJuti i%f» 

LD 

CTt1 

MU 

wv 

MU 

Goa 

MU 


DV 

Acc2 

MU 

cTi; 

AD 

TE 

MU 

Ceil U„ 

I 

ST 

(TB+DIP) 

ST 

KW 

MU 

CT12 

iSD 

lice.? 

ST 

DRIFT 

I'.lli 


LB LD 


AD 

H,, 
i r> 

ST 

\o 


Contd.. 

(continuod in the noxt ooltmn) 
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kl) 

BT 

SB 

t „ 
zo 

U'l’ 

(t.. 

‘ TP 

RS 

:.in 

fi 

CS 

• * 

AJi 

R 

P 

RS SB 

CT15 

3T 

Lgg2 

TP 

LB 

;.‘in 

kv.v,? 

LD 

E 

sp 

uv 


EN 

COS 

i4J 

H V/ 

5' 

ST 

Cos E 

sp 


(t.+Bf) 

IB 

Sin C 

37 

U Gin B 

•/> 

MU 

Sin S 

hitl 

li Gin n. 

ST 

Acc2 


nGv'/ 

IB 

Cos C 

Ei; 

GHHT 

MU 

Cos S 

S'r 

J.CU2 

SB 

Acc2 

Uj 

m 

MU 

Sin B 

Mli 

"t 

MU 

Cos E^j, 

LD 

‘'b 

Lvj'2 

ST 

Acc2 

3'i* 

li 

f’Z 

BN 

ARSIN 

',411 

CTf 

ST 

®sf 

Ki; 

ijism 

IB 

R^ 

i;iii 

Jfiii 

i,!i! 

ST 

CTB 

sm 

0T9 

MU 

ST 

MU 

MU 

^ 

WC 

RWC 

Cos 

CT14 


(Ounfcinuia in tho next colwnn) 


Gontd 


72 


ivHi 

Sin Sj, 

Bf 

E 



q. 

ST 


LD 

Sin s 

1.1) 

(TK+DIP) 

MU 

Cos C 

f.Ii’ 

Cp.'i 

ST 

Acc2 


ikCO . 

LD 

Cos 3 

ST 

nDv 

MU 

Sin C 

i.:n 

HWC 

iJ) 

Acc2 


Sin iHVf-’: 

MU 

Sin D 


CTIS 

DV 

Cos Eg^ 

4'fc i * 

I)i( 1 1'T 

EN 

ARSM 

ST 

<dil 

AD 

S 



P 

iiUl 

C ini S 

ST 

Aoc2 

ST 

An. cl? 

LD 

cDl' 

UJ 

l!l>V 

DV 

Cos 

r.tij 

Sin S 

AD 

Acc2 


A'fcC.HI 4f* 

ST 

T 

g . 

ST 

cOl * 

OU 

OUTPUT 1 

li) 

oDI 

LD 

E 

i.fU 

Sin S 

OU 

OUTPUT 2 

ST 

A«(32 

LD 

(t +dt) 

, z 

LD 

oOv 

OU 

OUTPUT 5 

Mil 

Cc«i 3 

LD 

R 

gz 

3h 

Aoc2 

OU 

OUTPUT 4 


®8f 

LD 

% 


(Continwd in the next colwm) 
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OU OOTPUT 5 

LD S, 

d 

OU OUTPUT 6 

Total N-umber of Instructions in main program = 276 


SUBPiODTIUE POE "SH" &"C0S" 


Entry for 

COS ST 

Acc2 

PB 

ST 

Acc2 



ID 

CT1 6 


LD 

CT 19 



SB 

Acc2 


ST 

SG 


Entry for 

SIW TP 

PA 

PD 

ID 

CT16 



iD 

CTI 7 


SB 

Acc2 



Ek ST 

Acc2 


TP 

PC 



ID 

CT16 


ID 

CT18 



SB 

Acc2 


SB 

Aoc2 



TP 

PC 


ST 

Aoc2 



ID 

Acc2 

PC 

ID 

Acc2 



SB , 

CT18 


ST 

THETL 



TP 

PB 


M 

Aoo2 



CS 

• * « 


CS 

m 0- • 



ST 

Acc2 


ST 

Acc2 



ID 

CT10 


DV 

CT20 



ST 

SG , 


AD 

CT10 , 



TP 

H) 


m 

Acc2 



(continued 

in the next 

column) 



Contd. 
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DV 

CT21 

AD 

CT10 

MU 

Acc2 

DV 

CT22 

AD 

CT10 

MU 

Acc2 

DV 

CT 25 

AD 

CT10 

MU 

THETA 

MU 

SG 

EM 

• • • 


Total number of instructions in Subroutine SBr/cOS 


SUBROUTINE SQHT . 
Algorithm used 8 



where Y = SQRT (c) 


ST 

Acc2 

ST 

AEG 

DD 

Acc2 

ST 

X 

0 

MU 

CT8 

ST 

Acc2 

ED 

AEG 


45 
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MU 

CT8 


ST 

Acc2 


LD 

AEG 


MU 

CT8 


m 

X 

0 



Acc2 


ST 

Acc2 


MU 

CT24 


ST 

EX 


ID 

X 

0 


SB 

ACC2 


TP 

PX 


CS 

4 « « 

EX 

SB 

EX 


TP 

IS 


ID 

Acc2 


EH 

• • « 

Total 

Ho. of ins 

trnctions 

SIB ROUTBIE 

4ESIH 



TP 

XA 


CS 

» « • 


ST 

AEG 


ID 

CT19 


ST 

SG 
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us 

• 0 • 

TP 

XC 


TP 

XB 

GS 

» » m 

lA 

ST 

km 

XC SB 

CT25 


ID 

CT10 

TP 

LP 


ST 

SG 

LD 

Acc2 

XB 

LD 

km 

MU 

SG 


MU 

km 

EU 

• • • 


ST 

X 





0 



LP 

ID 

EX 

Total No. 

of Ins true- 


MU' 

EX 




MU 

^0 




ST 

koc'Z 




LD 

EX 




IJ) 

CT10 




ST 

SG 




IiD 

CT10 




ST 

EX 




MU 

X1 




ST 

X1 




LD 

Acc2 




DV 

.XI 




ST 

Lee 2 




iiD 

ARG 




ST 

AEG 




LD 

Acc2 




(continued in the next column) 
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SYMBOL UBIE 


No. 

Symbol 

No. 

Symbol 

No. 

Symbol 

1 

DT 

2 

I* 

R 

3 

S 

i 

7 

\ 

5 

R 

P 

6 

s 

p 

T 

s 

8 

s 

9 

DID 

10 

E 

— S-'D. — ■ 

11 

B 

m: 

12 

w 

13 

E ' 

— a_ 

14 

s 

15 

dt 

16 

s 

17 


18 

Sin S 

P 

19 

Cos S 

]2_ 

20 

Sin S 

21 

Cos S 

22 

B 

Id 

25 1 


24 

W 

25 

C 

26 ' 


27 

W 

28 I 

Sin C j 

29~ 

Cos C 

50 

R w 

P 

31 


32 1 

Cos S^ 

33 

Sin S^ 

34 



wc 

35 


36 

: ■ "Vo 



1 

38 

R„ Cos D 

1 - - 

39 

Sin B 

i 

40 

D 

41 

T 

g 


Bwg 

43 

Sin Bwg 

44 

Cos Bwg 

1 - 

45 , 

w 

46 

TE 

47 

CT1 

48 

i GT2 

r 

49 

CT5 

50 

GT4 

51 

GT5 

_52 J 

CT6 

53 

CT7 

54 

GTS 

55 

CT9 

56 

CT10 

57 

CT11 

58 

CT12 

59 

CTI 5 

60 

GTI 4 

61 

0115 

62 

CT16 

_ . _ __ 

65 

GTI 7 

i 

1 

1 

1 

1 

1 

CT18 

0 

; 65 

j CTI 9 

i 

66 

CT20 


Contd. . . . . 
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SYMBOL TilELB (Contd.) 


No . 

Symbol 

No. 



Symbol 

Fo. 

Symbol 

67 

CT21 

68 

CT22 

69 

CT25 

70 

CT24 

71 

CT25 

72 

MliT 


















APISNPIX III 


CO NTROL circuits TO E XEGUIEE MlCEO/mCRO-OPSRi.TIOIIS 


A gate level design of the control circuits to execute the 
micro and macro-operations are given helow. The following notations 
are used to indicate various logical blocks. 


1) PP J 

2) L t 

3) D t 

4 ) A ! 

5) , E t 

6) C t 

7) CLKj 

8) MEXs 

9) RS s 

10) S t 

11) Lett 

12) MC I 
15) DC { 


Plip flop 
Latch 

Delay circuits 
AND gates 
OR gates 
Counter 
Clock pulse 
Multiplexer 

Reset flip flop/clear the register 
Set flip flop 
Addition complete 
Multiplication complete 
Division Complete 
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M1 ! 


o 8 



0 8 

When, the M1 line is pulsed, latch I4 clears IviAR. The input 


pulse is delayed by D3 till the MA.R is cleared. Output of DJ is given to 
Mh gates A3O through A38 vdiich transfer the contents of LC to B/iiR. 



0 8 

?Jhen the input line of M2 is pulsed, D4 delays the pulse till 


L5 clears M/iR. Output of D4 enables transferring of the contents of the 
address part of IR to MR by AND gates A39 through M7* 
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;i 'Toggle^ all flip flops 

.* ii-CC **1 

__ L_i 

then the input line of Wf is pulsed, L10 pulses the J and K 
inputs of all flip flops of Acc ^diich toggles all flip flops in Acc , 
effecting 1'® complementation. 
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goggle ' all flip flops 
~ 1 __ 


o 


47 



I’he-n the input line of M 9 is pulsed, FP6 is set and the 'n' 
counter is reset. 1 output of FP6 connects i-cc and B registers to the 
Adder through ANT) gates AI 58 thro-ugh A 5 OI and the Addition commencei^. 
The 'n ' coxxnter C2 starts counting the pulses. Ihen the counter has 
counted upto 'n', output of 02 resets PP6 and gives addition complete 
(AC) signal to the MSG. 
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Wien input line of M16 is pulsed, FF10 is set and the 'n’ cotmtor 
is rest. 1 output of FF10 connects the K and bits 4 throxjgh 11 of IH to the 
9 bit parallel Adder, through the AHl gates A861 through A877 , addition 
commences and C5 starts counting- the clock pulses. Fhen the counter has 
counted upto ’n’, output of C5 resets FF10 and gives addition complete 
signsil to MSG8. Here LC has 9 bits v^here as the signed displacement i,a. 
bits 4 through 11 of IR is of 8 bits. Hence bit 4 of IH is connected to 
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MixO 'll 


Inputs 



When the input line of MCI is pulsed, FPU is set and C6 is 
nesot. 1 output of PP11 enables the 1 out of 16 selector to select the input 
line as indicated by bits 8 thro-ugh 11 of IR and connects it tc the A/D 
converter. The C6 starts counting. ?hen C6 has counted uptc 'nS where 
'n' is the number of clock time A/d cnversicn takes, output of 06 resets 
PP11 and shifts the contents of l/O Buffer Re^gister to Iicc. 
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ISTien input line of IvI/iC2 is pulsed FF12 is set, Of is reset and the 
eontents of i-cc are shifted to l/O Buffer register-. 1 output of PP12 enables 
the B/A converter to convert the contents of l/O Buffer Register into ii-nalog 
signal. The 8 line Multiplexer (MK) connects the output of the D/A converter 
to the output line indicated by bits 9 through i1 of IR. '^hen the D/A 
•onversion is over, output of C 7 resets FP12 and gives end of operation signal. 
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GLOSSARY 


Acc t Accijmulator Register 

B^s Future bearing of target 

B^ I Rate of change of bearing in lateral plane 

B^^^J Relative Wind bearing 

C^s Combined drag coefficient 

LB: Error in prediction of bearing 

LEs Error in prediction of elevation 

DR: Error in prediction of range 

B^: Predicted future elevation 

E : Gun elevation 

g 

F^: Improvement factor for bearing. 

: Improvement factor for elevation 

Improvement factor for range 
IC i Integrated circuit 

IR: instruction Register 

LC: Location Counter 

LSI: Large Scale integration 

MAE: Ifemory Address Register 

ICBR: Ifemory Buffer Register 

MQ: Multiplier Quotient Register 

MSG: Micro-Sequence Generators 
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Yi 


Mi 


Mediim Scale Integration. 

Mean time "between failnres 

Brojective velocity 

Itedicted future range 

Present range 

Pate of change of range 

Rotation, of sight line 

Present elevation 

Rate of change of elevation 

Small Scale Integration 

G-un training 

Present bearing 

Transistor-transistor logic 

"Wind velocity 

ingle of approach 

Diving angle of aircraft 

Sampling to interval 

Acceleration due to gravity 

Predicted time of flight of the Shell 

of to-rg-et- 

Angular rate of the target 



